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The scaled factorial moments in second-order quark-hadron phase transition are reexamined 
within the Ginzburg-Landau description, with the spatial fluctuations of phase angle of the complex 
held <p taken into account rigorously. Scaling behaviors between F q and F 2 are shown, and the 
exponent v is found very close to the one without spatial fluctuations. 
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The idea of scaled factorial moments [1] has been used since being proposed to study dynamical fluctuations 
extensively in both theoretical investigations and experimental analyses due to the powerful ability to eliminate 
statistical (non-dynamical) fluctuations. Some scaling laws for the moments have been predicted theoretically and 
observed in high energy hadronic and nuclear interactions [2], In high energy heavy-ion collisions, a new matter state 
-quark-gluon plasma (QGP)- may be produced. The newly produced hot system will cool with its expanding and 
will undergo a phase transition from the deconfined QGP to confined hadrons which will be detected in experiments. 
Since the existence of the phase transition is associated with properties of the nontrivial chromodynamical vacuum, 
the study of quark-hadron phase transition has been a hot point in both particle physics and nuclear physics for 
more than a decade. Besides the unique features of quantum chromodynamics, the lack of control of the temperature 
in the phase transition distinguishes the problem from the standard critical phenomena such as ferromagnetism. 
The nonperturbative nature of hadronization process in the phase transition precludes at this stage any observable 
hadronic predictions from first principles. Hence, among many phenomenological models, the Ginzburg-Landau 
description has been used as a framework to make predictions, which can be compared with experiment [3, 4], on 
the hadronic observables. Up to now, this model has been used in the study of various scaling behaviors in the 
phase transition, such as those of the scaled factorial moments in both first-order [5-7] and second-order [8,9] phase 
transitions, the multiplicity difference correlators [10-12], the multiplicity distributions [13,14], etc. 

In the Ginzburg-Landau description of second-order phase transition, the scaled factorial moments can be expressed 
as [9] 

Fq = fq/fl , fq = ^ j Ff) dz \ (j) [^ eX P (-F[0]) (1) 

with Z = J T><j)exp(—F[(j)}), the free energy functional F[(j)\ = f dz[a \ <j) | 2 +b \ <f> | 4 +c | \j<f | 2 ], a oc (T — Tc) 
representing the distance from the critical point, b and c larger than zero. Here \<f>\ 2 is associated with the multiplicity 
density of the system. Similar expressions can be derived for other quantities mentioned above. In all former studies of 
second-order phase transition the gradient term is simply taken to be zero, i.e., the field <f> is regarded spatially uniform. 
The spatial integral of the functional over a two-dimensional bin with size 5 2 is then F[(j>] = <5 2 (a, | <^> 2 | +b \ <fi | 4 ) 
(Calculations based on this functional will be referred to mode 1 in this paper). This is of course a very crude 
approximation. The advantage of such an approximation is, however, that it turns the functional integration into a 
normal one. Thus, the calculation becomes quite easy under the approximation. Numerical results do not show the 
so-called intermittency behavior, but the F-scaling, F q oc F^ with universal scaling law (3 q = (q — l) 1 ', is shown to 
be valid. The exponent v is called as a universal one in the sense that it is insensitive to the values of the parameters 
in the model and that it is completely determined by the structure of the functional concerned. 

The contribution from the gradient term to the moments and the exponent v should be evaluated in some way. Once 
the gradient term is taken into the functional, one is faced with serious difficulty in the calculations, considering the fact 
that the value of parameter b for the </> 4 term can be determined in no way from first principles or from experimental 
input and may be very large. Even the parameter b is indeed very small, negative value of a in our interested region 
also excludes the possibility of performing usual perturbative calculations. The role played by the gradient term is 
investigated in Refs. [15, 16]. In Ref. [15] (f> in each bin is still uniform, but the values of cj) in all neighboring bins 
are taken to be </>o, held configuration corresponding to the minimum of “potential” V(<j>) = a \ 4> \ 2 +b \ <f> | 4 . So 
the square of the gradient of <f> is 5~ 2 (<f> — <f>o) 2 - This approximation also transforms the functional integration into a 
normal one. Numerical results show that the universal scaling law /3 q = (q — 1)" is still valid and that the exponent 
v is almost the same as without the gradient term. In Ref. [16] the details of spatial fluctuations of in a bin is 
simulated by the Ising model for one-component spins s. Each bin is assumed large enough to contain several spin 
sites. The multiplicity in a bin is associated with the sum of spins on all sites in the bin by 
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Hi = A I J2 s i I 2 0 ( Y1 s j) > 

j’Gbin j'Gbin 

in which scale factor A relates the hadron density to the lattice spin density in the Ising model. Since A relates a 
physical space to a mathematical space and there is no way to determine it a priori, several A’s are used in Ref. [16]. 
This time, the exponent v depends on the unknown temperature, and, after averaging over the temperature, v is still 
in the range given in Refs. [8, 9, 15]. 

Though the simulation in Ref. [16] is convincing, it is for lattice with one-component spins. In the Ginzburg- 
Landau model for second-order phase transition, the field </> is a complex number, or in other words, (f> has two 
components. Thus at first glirnps the simulation in Ref. [16] does not correspond to the real problem discussed in 
the Ginzburg-Landau model. But as will be shown soon the simulation relates to the physics in an indirect way. The 
simulation, however, is still useful because one can translate the two-component field into a one-component one, as 
will be illustrated below in this paper. 

In this paper, it is tempted to investigate the universality of the exponent v 1 with the spatial fluctuations of the 
phase angle of the complex field </> fully taken into account. The magnitude fluctuations of </> is neglected in this 
investigation for reasons which will be given below. As will be seen soon, the contribution from spatial fluctuations 
of the phase angle of the field <f> can be taken into account in a complete and rigorous way, and the integration over 
the spatial fluctuations of the phase angle of the field <j> will reduce the problem to one with one-component field. 

Our first observation is that all terms except the gradient term in the functional integral of Eq. (1) depend only 
on | 4> | 2 . Then it is convenient to write the two-component field as a complex number in the form <j> = </>r exp(?’(/i). 
The spatial fluctuations of the field can be those of the magnitude <j> R and/or of the phase angle <f >i (or orientation in 
an abstract space). The gradient term turns out to be 

| v</> | 2 = (V</r) 2 + </r(V</>i) 2 • (2) 


If v0i is fixed in the process (either because of the uniformity of (j>\ or uniform change of (j> i), the effect of (\/(f> i) 2 can 
be effectively considered by a new parameter a 1 in place of a with a 1 = a + c(v<^i) 2 , reducing the problem to former 
case if the \/4>r term is also neglected. Needless to say, such a uniformity can be realized only in the case with very 
strong correlation between particles within the bin and, most probably, does not correspond to our case under study. 
Generally, the phase angle <fii can be in any form, and the full contribution due to its fluctuations must be evaluated. 
Fortunately, the integral over <f>\ can be carried out easily since it is of Gaussian form. Using a formulus 

J 4>RD(j)RD(j)\ exp[— J dx(jy^{\/4>i) 2 ] = const • j Vcj) r , 

one transforms the two-fold functional integral into a one-fold one. The unimportant constant will be cancelled in the 
expression for f q . Then one has 


, = I V( t>K (Is dz^) 9 exp(-F[(/> R \) 

/ exp(—E[</ r]) 

with functional E[</r] exactly the same form as the original F [</>]. The important difference between this expression 
from Eq. (1) is that the functional integral variable in this new expression is a real function instead of a complex 
function in Eq. (1). Similar expression can also be derived for the case of first-order phase transition, with an extended 
functional F[</>r]. 

Now we take the field </r (magnitude of (f >) uniform, or in other words, the gradient term of </>r is omitted. 
(Calculations based on this approximation will be referred to mode 2 in this paper.) Based on the work Ref. [16] one 
can drop off the V</r term, because the problem now is exactly within the scope of Ref. [16], and the conclusions 
in Ref. [16] encourage us to neglect the spatial fluctuations of </r as long as the universal scaling exponent v is 
concerned. Then one gets the factorial moments as functions of variable x 

, = / 0 °° dyy 2q exp (xy 2 - y 4 ) 

' q f 0 °° dy exp (xy 2 - y 4 ) 

with x = ad 3 / 2 /^ 1 / 4 . From this expression the scaled factorial moments In F q can be calculated, and the results are 
shown as functions of — In a; in Fig. 1 for q from 2 to 8 within the range x £ (0.5,4.0). One can see clearly that no 
strict intermittency can be claimed since all F q approach finite values in the limit x —> 0. So, no intermittency is 
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shown in the phase transition, as shown in former studies. More importantly, the power law can be found between 
F q and A 2 , as shown in Fig. 2 with the same data as in Fig. 1. 

For the convenience of comparison with former case, we write down the expressions of the scaled factorial moments 
without spatial fluctuations (mode 1 in this paper), which can be read 

, _ /o°° dyy q exp (xy - y 1 2 ) 
fo° dyexp(xy — y 2 4 5 6 7 * ) 


with x = aS/Vb. Numerical results for In F q in this mode are shown in Fig. 3. In the upper part of the figure In F q 
are shown as functions of — In 2 : for q from 2 to 8 with x in the same interval x € (0.5,4.0), and in the lower part 
InFq are shown as functions of In A 2 with the same data as in upper part. One can see from upper part of the figure 
that the general behaviors of In F q as functions of — In x is similar to those in Fig. 1, though the definition of x in 
this case is different from that for Fig. 1. The values of In F q in the two cases are also different. For same value of x, 
In F q in the former case have larger values. This difference is reasonable if one notices the difference in the definition 
of variable x. What interests us is the scaling law between F q and F 2 . The power law scaling between F q and F 2 can 
be seen obviously in the lower part of Fig. 3, the same as shown in other studies cited in the references. 

From Fig. 2 and the lower part of Fig. 3, one can get the scaling exponents f3 q for the two different modes by fitting 
the curves. (3 q can also be given analytically. One can expand the expressions for In F q in the two modes as power 
series of x in small x limit, and then one gets the slopes K q for lni^ and (3 q = A' g /A" 2 . If the scaling is good enough 
as shown in this paper, the exponents obtained from these two different ways can be regarded equal. The expressions 
for K q for the two modes in this paper are 


K q 

K<1 


r(g/2 + i) £(3/2) r(i) 

r ( g /2 + i/ 2 ) q r(i) [q j r(i/ 2 ) 

F(g/2 + 3/4) r(5/4) , r(3/4) 

r(g/2 + l/4) 9 r(3/4) [q j r(l/4) 


for mode 1 , 


for mode 2 . 


One can find only a small difference between the exponents v from these two expressions. The results are shown in 
Fig. 4. In mode 1 (without spatial fluctuations) ^=1.3335, and in mode 2 (with spatial fluctuations of the phase angle 
of the held </>) i/=1.2772. The exponents obtained from these analytical expressions are very close to the ones from the 
fitting. The universal exponent v are also very close to one another and can be regarded as the same within accuracy 
4%. Physically, these two modes correspond to different situations. In mode 1 no spatial fluctuation of 4> is in the 
problem, but in mode 2 the spatial fluctuations of the phase angle of the complex held 4> are fully evaluated. Since 
these two different considerations give very close exponents v, one can say that the exponent v is indeed insensitive 
to the spatial fluctuations. It should be pointed out that the slight difference of the exponent v for mode 1 in this 
paper from that in Ref. [9] comes from the different x regions from which /3 q is calculated. In Ref. [9] f3 q (thus v) is 
obtained by fitting in larger x, but in this paper (3 q corresponds to small x region. 

In summary, the scaled factorial moments in second-order quark-hadron phase transition are reexamined with the 
Ginzburg-Landau description, with the spatial fluctuations of the phase angle of the held fully taken into account. 
Scaling behaviors between F q and A 2 are shown, and the exponent v is found universal within accuracy 4%. It will 
be interesting to investigate the effect of the full spatial fluctuations of the held on the exponent v. This topic will be 
left for further discussions. 
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Figure Captions 

Fig. 1 Dependences of In F q on the bin width — In x after the contribution from spatial fluctuations of the phase angle 
of the field fully taken into account (mode 2). Curves from lower to upper are for q from 2 to 8, respectively. 

Fig. 2 Scaling behaviors of In F q vs lni 7 ^ for the same data as in Fig.l. 

Fig. 3 Upper part: In F q as functions of — In a: without spatial fluctuations (mode 1); Lower part: Scaling behaviors 
between \nF q and lni 7 ^, with the same data as in the upper part. 

Fig. 4 Scaling behaviors of In /3 q as function of ln (<7 — 1) for the two modes. 
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